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A common signaling method used in the 
telephone network is dual-tone multifrequency 
(DTMF) signaling. In this scheme, pairs of tones 
are used to signal the digits through 9, 
pound (#), star (*), and the digits A, B, C, and 
D. For each pair, one of the tones is selected 
from a low group of four frequencies, and the 
other from a high group of four frequencies. 
The correct detection of a digit requires both a 
valid tone pair and the correct timing intervals. 

DTMF signaling is used both to set up a call 
and to control features such as call forwarding 
and conference calling. In some applications, it 
is necessary to detect DTMF signaling in the 
presence of speech, so it is important that the 
speech waveform is not interpreted as a valid 
signaling tone. Standard test procedures have 
been published [1]* to verify the performance of 
DTMF receivers for valid tones and for the 
recognition of speech signals as tones. 

The implementation of a DTMF receiver involves 
the detection of each of the signaling tones, 
validation of a correct tone pair, and timing to 
determine that a digit is present for the correct 
amount of time and with the correct spacing 
between tones. In addition, depending on the 
algorithm used to detect frequencies, it is 
sometimes necessary to perform additional 
tests to improve the performance of the 
decoder in the presence of speech. Current 
DSP technology allows several DTMF receivers 
to be implemented on a single device. A DSP 



[ ] indicates a reference listed at the conclusion of this 
application note. 
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implementation is useful in applications in 
which the digitized signal is available and 
several channels need to be processed, such 
as in a private branch exchange (PBX), or 
where other functions can be included on the 
same device. 

This application note describes a DTMF 
receiver based on the discrete Fourier 
transform (DFT). Using this algorithm, 16 DTMF 
receivers can be implemented on a WE DSP32 
Digital Signal Processor using internal RAM (for 
both data and program memory) and executing 
at a 160 ns instruction cycle time. The input 
data is assumed to be //-law coded, and is 
received on the DSP32 serial interface. 
Detected digits can be read from the DSP32 
parallel port. 

DTMF Requirements 

Figure 1 shows the matrix of frequencies used 
to encode the 16 DTMF symbols. Each symbol 
is represented by the sum of the two 
frequencies that intersect the digit. The row 
frequencies are in a low band, below 1 kHz, 
and the column frequencies are in a high band, 
between 1 kHz and 2 kHz. The digits are 
displayed as they would appear on a 
telephone's 4x4 matrix keypad (on standard 
telephone sets, the fourth column is omitted). 
DTMF receivers are required to detect 
frequencies with a tolerance of ±1.5 percent as 
valid tones. Tones that are offset by ±3.5 
percent, or greater, must not be detected. This 
requirement not to detect tones is necessary to 
inhibit the detector from falsely detecting 
speech and other signals as valid DTMF digits. 
The receiver is required to work with a worst- 
case signal-to-noise ratio (SNR) of 15 dB, and 
with an attenuation of 26 dB. 

Another requirement is the ability to detect 
DTMF signals when the two tones are received 
at different levels. The high-band tone may be 
received at a lower level than the low-band 
tone due to the attenuation characteristics of 
the telephone network. This level difference is 
called twist, and the situation described is 
called normal twist. Reverse twist occurs when 
the low-band tone is received at a lower level 
than the high-band tone. The receiver must 
operate with a maximum of 8 dB normal twist 
and 4 dB reverse twist. 

In addition to the frequency, noise, and twist 
requirements, the DTMF signal must meet 
timing requirements for duration and spacing 



of digit tones. Digits are required to be 
transmitted at a rate of less than ten per 
second. A minimum spacing of 50 ms between 
tones is required, and the tones must be 
present for a minimum of 40 ms. Any tone- 
detection scheme used to implement a DTMF 
receiver must have a significant time resolution 
to verify correct digit timing. 

A final requirement for the receiver is that it 
operates in the presence of speech without 
incorrectly identifying the speech signal as a 
valid DTMF symbol. This is referred to as talk- 
off performance. Although this requirement is 
not stated in strict numerical terms, standard 
recordings such as the Mitel DTMF test tape [2] 
contain speech segments that are used to test 
the receiver's performance under these 
conditions. 

The remainder of this application note 
describes the algorithm, implementation, and 
performance of the DFT-based WE DSP32 
Digital Signal Processor DTMF receiver. 

1 209 Hz 1 336 Hz 1 477 Hz 1 633 Hz 



697 Hz 



770 Hz 



852 Hz 



941 Hz 



Figure 1. 4x4 Matrix Telephone Keypad 

The Basic Algorithm 

The general approach taken by this algorithm 
for DTMF tone detection is to take the Fourier 
transform of the observed signal and search 
for energy at the frequencies of interest. Since 
the algorithm is implemented on a DSP, a 
Discrete Fourier Transform (DFT) is used. The 
analysis frame must be long enough to resolve 
the DTMF frequencies, but short enough to 
detect the minimum length tone. A 12.75 ms 
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frame at a sampling rate of 8 kHz is a good 
choice [3]. 

To guard against talk-off, the total energy is 
calculated and compared to the sum of the 
largest magnitudes in the high- and low- 
frequency bands. If a valid tone pair is being 
transmitted, then the two values should be 
equivalent. But if there is speech present, the 
total energy should be much greater. 

The algorithm also exploits the rich harmonic 
structure of the speech signal to assist in 
guarding against talk-off. If speech has energy 
at a frequency, f-i , it most likely has significant 
energy at twice the frequency, 2^ . Simulation 
of a DTMF tone by speech is prevented by 
calculating the energy at each tone and at its 
second harmonic. If the second harmonic is 
below a threshold proportional to the 
fundamental, a pure tone is present. Otherwise, 
it is assumed that a speech signal is present. 

In calculating the DFT, the Goertzel algorithm 
[4], a method for calculating any single 
coefficient of a DFT, is chosen over a fast 
Fourier transform (FFT) algorithm. There are 
two reasons for this. In order to obtain the 
required frequency resolution at an 8 kHz 
sampling rate, a 256-point FFT would be 
required. Since the algorithm for tone detection 
requires knowledge of the energy at only 16 
frequencies, it is more efficient to execute the 
Goertzel algorithm for these 16 frequencies 
than to perform the FFT for all 256 frequencies. 
In addition, the Goertzel algorithm is recursive, 
eliminating the need to store 256 samples for 
the FFT for each DTMF detector. This saves 
both real-time and data memory. 

The Goertzel algorithm can be thought of as a 
matched filter for each DFT coefficient. The 
transfer function for the filter is: 



H k (Z)= 



1-WfcZ" 1 



1-2cos(-^)Z- 1 +Z" 2 



N 

where W N = exp(-j27r/N), and N is the length of 
the observation window for the DFT. The flow 
graph of this transfer function is shown in 
Figure 2. 

Initially, the state variables of the filter are set 
to zero. Then, the filter is executed N times. The 
output at this point, y k (N), is the k'th coefficient 
of a length N DFT. Notice that the filter is 
implemented as a direct form II second-order 



section. The recursive part of the filter is on the 
left-hand side of the delay elements, and the 
nonrecursive part is on the right. Since only 
the output at time N is needed, it is only 
necessary to compute the nonrecursive part of 
the filter after the last iteration of the recursive 
part. A further simplification in the algorithm is 
made by observing that only the square of the 
magnitude of the DFT coefficient is needed. The 
nonrecursive calculation of the DFT coefficient 
is: 



y k (N)=S(N)- 



-W^S(N- 



D 



where S(N) and S(N-1) represent the value of 
the state variables at times N and N-1. It can 
be shown that: 



|y k (N)| ; 



,27rk 



| S(N) | 2 -2cos(-^)S(N)S(N-1)+ 1 S(N-1) | 2 

Therefore, it is only necessary to store the 
value, 2 cos (27rk/N), for each coefficient to be 
evaluated. 




S(n-2) 



Figure 2. Flow Graph of the Transfer Function 

Implementation 

The flow chart for the DTMF tone detector is 
shown in Figure 3. 

At the beginning of each frame, the state 
variables of each of the 16 Goertzel filters are 
set to zero. Then, for 12.75 ms (102 samples at 
a sampling frequency of 8 kHz), the recursive 
part of the Goertzel algorithm is executed. At 
the end of each frame, the square of the 
magnitude of the coefficients at each frequency 
are calculated. The following five tests are then 
performed on these 16 values: 

1. The Magnitude Test In each group of four 
tones (the row and column tones), the 
frequency whose coefficient has the largest 
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magnitude is found- This frequency is 
tentatively referred to as the detected 
"tone" for that band. The detected tones' 
magnitudes are compared to a threshold. 
If either tone's magnitude is smaller than 
the threshold, then this test has failed. 

2. The Twist Test. Twist is the ratio of the 
energy in the lower-band tone to the 
energy in the upper-band tone. If the 
measured twist falls outside the region of 
-4 dB to +8 dB, the test has failed. 

3. The Frequency Offset Test. The energy of 
the largest tone in each band is compared 
to the energies of the other frequencies in 
that band. If the difference falls below a 
threshold in either band, the test fails. 

4. The Tone-t©-TotaI Energy Test. Let d , c2, 
and c3 be three different constants, each 
greater than one. The energy of the low- 
band detected tone is weighted by d , the 
energy of the high-band detected tone is 
weighted by c2, and the sum of the low- 
and high-band tones' energies are 
weighted by c3. If any one of these terms is 
smaller than the total energy, then the test 
has failed. 

5. The Harmonic Rati© Test. The energy of 
the detected tone is compared to the 
energy of its second harmonic. If the 
energy of the second harmonic is not 
sufficiently smaller than the energy of the 
detected tone, the test has failed. 

If all of the above tests pass, the tone pair is 
decoded as an integer between and 15, 
representing the digits through 9, A through 
D, * and #, respectively. This value is placed in 
a memory location designated D(k) and is the 
digit at frame k. If any of the tests fail, then -1, 
representing "no detection," is placed in D(k). 
For a new valid digit to be declared and sent to 
the parallel port, D(k) must be the same for two 
successive frames. If it is valid for more than 
two successive frames, the receiver is detecting 
the continuation of the previously validated 
digit, and a new digit is not output. 

Performance 

The DTMF detector described above has been 
verified by using the Mitel DTMF test tape. 

This section documents the results for each 
test included on the tape. 

Test 1 - This test determines the frequency 



offset range of the receiver. Acceptable 
performance requires detection of the tone 
when the magnitude of the frequency offset is 
less than ±1.5 percent. When the magnitude of 
the frequency offset is greater than ±3.5 
percent, the receiver is required not to detect 
the tone (see Table 1). 

Table 1. Frequency Offset Range of the 
Receiver 



Frequency 
Deviation 


Percent Deviation at Cutoff 

Digit 

1 5 9 D 


Low group 


2.4 2.4 2.7 2.7 
-3.5 -2.7 -2.2 -2.3 


High group 


2.3 2.1 2.0 1.8 
-2.4 -2.0 -1.9 -1.7 



Test 2 - This test verifies that the receiver 
operates with a normal twist of up to 8 dB and 
a reverse twist of up to 4 dB (see Table 2). 

Table 2. Verification of Receiver's Twist 
Operations 



Twist 


Attenuation at Cutoff (in dB) 

Digit 

1 5 9 D 


Reverse 


5.3 4.3 4.4 4.1 


Normal 


8.2 8.1 8.3 8.4 



Test 3 - Dynamic range test. A receiver must 
detect touch-tones with a dynamic range of 26 
dB. The receiver detects touch-tones with a 
dynamic range of 34 dB. 

Test 4 - Guard time test. A receiver must 
detect DTMF pulses as short as 40 ms. The 
receiver detects pulses as short as 34 ms. 

Test 5 - Signal-to-noise ratio test. A receiver 
must detect touch-tones with a SNR of 15 dB. 
The receiver detects all touch-tones at this 
level. 

Test 6 - Talk-off test. Acceptable talk-off 
performance requires less than 30 false digit 
detections in the talk-off section of the Mitel 
tape. The receiver detected 5 false digits. 

Hardware Description 

The WE DSP32 Digital Signal Processor 
receives its sampled data input from a time- 
division multiplexed (TDM) serial bit stream 
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which contains 16 channels of 8-bit //-law 
samples. A sample timing diagram is shown in 
Figure 4. The rate of the input load clock (ILD) 
must be 128 kHz (one clock cycle for each 
channel). The input clock (ICK), which is the 
rate at which the bits are input to the DSP32, 
should equal 2.048 MHz. The signal SY, which 
is the external synch pulse, tells the DSP32 
when the first channel is being transmitted. 
When the ILD corresponding to channel 16 



goes high, the SY should go low coincidental^ 
and go high 800 ns before the ILD 
corresponding to channel 1 goes high. The 
DSP32 DMA hardware handles the loading of 
input data to a buffer in RAM. 

The output of the DSP32 is connected to an 
external microprocessor. When a digit has 
been validated, it is written to the lower 8 bits 
of the parallel I/O data register (PDR). The 
lower 4 bits contain the decoded digit, and the 




Figure 3. Flow Chart for the Dual-Tone Multlfrequency Detector 
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upper 4 bits of the byte contain the channel 
number. When the PDR is written, the Parallel 
data full (PDF) pin goes high, which can be 
used to interrupt the external microprocessor. 
The microprocessor then has a maximum of 40 
n$ to read the data before it is overwritten. 

Source Code 

Program 1 implements 16 DFT-based DTMF 
receivers on the AT&T WE DSP32 Digital Signal 



Processor operating with a 160 ns instruction 
rate. Both the program and data are stored in 
on-chip RAM. The state variables for the DFT 
are in bank 1, and everything else is in bank 0. 
During the 125 //s sample period, there is time 
to execute the 16 receivers and an additional 
130 instructions. 
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Figure 4. DSP32 TDM Data and Microprocessor Interfaces 
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Program 1. DTSVIF Receiver Source Code 



. global _iirstvr ,_iirincr ,_ttdigit , inbuf , buf io 



dauc=0 


ixO 




ioc=0x 


:40 




rl5=12 






rl6 = 4 






begin:if (syc) 


goto ] 


begin 


nop 






ioc=0x 


:2040 




pin=inbuf 




rl7=-60 




r5=102 


:-2 




r6=energy 




rlO=_iirzero 




r7=16- 


2 




clrnrgrif (r7- 


- >= 


) goto clrnrg 


*r6 + + 


= aO = 


*rlO 


r7=32- 


2 




r2=_iirco 




r3=_iirstvr 




clear :*r3++ = 


aO = *rlO 


*r3+ + 


= aO = 


*rlO 


*r3+ + 


= aO = 


*rlO 


*r3 + + 


= aO = 


*rlO 


*r3+ + 


= aO = 


*rlO 


*r3 + + 


= aO = 


*rlO 


*r3 + + 


= aO = 


*rlO 


*r3+ + 


= aO = 


*rlO 


*r3 + + 


= aO = 


*rlO 


*r3+ + 


= aO = 


*rlO 


*r3 + + 


= aO = 


*rlO 


*r3 + + 


= aO = 


*rlO 


*r3 + + 


= aO = 


*rlO 


*r3 + + 


= aO = 


*rlO 


*r3 + + 


= aO = 


*rlO 


if (r7 


r > = 


0) goto clear 


*r3 + + 


= aO = 


*rl0 



/* 8-bit input , passive elks, DMA on IBF */ 



/* wait for sync to go high before */ 
/* beginning DMA */ 

/* 8-bit input, passive elks, DMA on IBF */ 



/* DFT frame size - 2 */ 



/* 16 receivers - 2 */ 

/* zero the energy locations on a */ 

/* frame by frame basis. */ 

/* pointer to coefficients */ 

/* pointer to state variables */ 



/* zero the state variables on a */ 
/* frame by frame basis. */ 

/* */ 

/* This section sets up the DMA for each receiver. After 16 samples */ 

/* (one for each receiver) have been received, they are converted to */ 

/* floating point representation and stored at bufio. DMA is then */ 

/* re-initialized to continue receiving data. */ 

/* */ 

sample: pin-ilimit 

iwait:if (le) goto iwait /* wait for input data to be gathered */ 

pin-ilimit 

pin=inbuf 

rl=inbuf 

r4=buf io 

r6=energy 

r3=14 /* 16 receivers -2 */ 

ioxfr-.if (r3 — >= 0) goto ioxfr 

*r4++ = al = ic(*rl++rl6) /* move data from inbuf to bufio */ 

rll=14 /* 16 receivers - 2 */ 
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/* pointer to state variables */ 

/* stores updated state variables */ 



/* 2nd state var . d[l,l] */ 
/* 1st state var. d[l,2] */ 



/* */ 

/* This section computes the LHS of the DFT for the 1st and 2nd order */ 

/* harmonics of each tone for each sample. It is performed 16 times, */ 

/* once for each receiver. The state variables for each receiver are */ 

/* stored in a sequential stack on top of each other. */ 

/* */ 

rl=buf io 

r3=_iirstvr 

r4=r3+4 
dft:al=*rl 
/* 697hz 1st harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2++ 

*r4++rl5 = aO = aO - *r3++ 
/* 770hz 1st harmonic LHS */ 

aO = al -tr ( *r4 — = *r3 + +) * *r2 + + 

*r4++rl5 = aO = aO - *r3++ 
/* 852hz 1st harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2++ 

*r4++rl5 = aO = aO - *r3 + + 
/* 941hz 1st harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2 + + 

*r4++rl5 = aO = aO - *r3++ 
/* 1209hz 1st harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2 + + 

*r4++rl5 = aO = aO - *r3 + + 
/* 1336hz 1st harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2++ 

*r4++rl5 = aO = aO - *r3 + + 
/* 1477hz 1st harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2++ 

*r4++rl5 = aO = aO - *r3 + + 
/* 1633hz 1st harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2++ 

*r4++rl5 = aO = aO - *r3 + + 
/* 697hz 2nd harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2++ 

*r4++rl5 = aO = aO - *r3+ + 
/* 770hz 2nd harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2++ 

*r4++rl5 = aO = aO - *r3 + + 
/* 852hz 2nd harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2++ 

*r4++rl5 = aO = aO - *r3++ 
/* 941hz 2nd harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2 + + 

*r4++rl5 = aO = aO - *r3++ 
/* 1209hz 2nd harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2 + + 

*r4++rl5 = aO = aO - *r3++ 
/* 1336hz 2nd harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2 + + 

*r4++rl5 = aO = aO - *r3 + + 
/* 1477hz 2nd harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2++ 

*r4++rl5 = aO = aO - *r3+ + 
/* 1633hz 2nd harmonic LHS */ 

aO = al + (*r4 — = *r3++) * *r2++rl7 

*r4++rl5 = aO = aO - *r3+ + 

if (rll — >= 0) goto dft /* 16 receiver counter */ 

*r6++ = a2 = *r6 + al * *rl++ /* update energy */ 

if (r5-- >= 0) goto sample /* 205 sample frame counter */ 

nop 



/* d[i,i 

/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 

/* d[l,l 
/* d[l,2 



*/ 

V 

V 

*/ 

V 

■*/ 

*/ 

*/ 

V 
V 

V 

*/ 

*/ 

*/ 

*/ 
*/ 

*/ 
*/ 

*/ 
*/ 

*/ 
*/ 

*/ 
*/ 

*/ 

*/ 

*/ 
*/ 

*/ 
*/ 
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/* */ 

/* This section computes the RHS of the DFT for the 1st harmonics */ 

/* of each tone. The magnitude is squared and stored in consecutive */ 

/* locations starting at _iirstvr. */ 

/* */ 

ioc=0x40 

rl6 = 

rl2=energy 

rl8=_iirstvr 

rl9=_ttdigit 

rll=14 /* 16 receivers - 2 */ 

timing :r2=_iirco 

r4=rl8 

rl5 = 8 

rl = 6 

r3=rl8 /* points to location of mag~2 of each tone */ 

right:a2= *r4 * *r4 + + /* Sn * Sn-1 */ 

al=*r4 — * *r4 /* Sn-1~2 */ 

aO = al + *r4++rl5 * *r4 /* Sn~2 + Sn-1~2 */ 

if (rl — >= 0) goto right 

*r3++ = aO = aO - a2 * *r2++ /* |y|~2 */ 

/* V 

/* This section finds the largest value in the low and high frequency */ 

/* band and stores the 2 values at _iirstvr+32. It also computes their */ 

/* location in the low and high frequency band and stores these 2 */ 

/* values at iirincr. */ 

/* " */ 

r8=_iirincr 

r4=rl8+32 /* loc. of largest mag. in each band */ 

r3=rl8 /* loc. of mag~2 of each tone */ 

r7 = 
lrgmag3 : r6=l 
r5 = 
rl = 
r2 = 
*r4=a0=*r3++ /* writes possible large mag. temporarily *r4 */ 

lrgmag4 : al=*r3+ 

nop 

a2=al-a0 

3*nop 

if (ale) goto lrgmagl /* test which value is larger */ 

nop 

*r4=a0=al 

r2=r2+l 

rl=r2 

goto lrgmag2 

rl=rl+r5 /* contains the row & col increments */ 

lrgmagl :r5=r5+l 
lrgmag2:if (r6 — >= 0) goto lrgmag4 

nop 

* r 8++=rl /* contains the row & col increments */ 

if (r7 — >= 0) goto lrgmag3 
r4=r4+4 
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/* */ 

/* This section compares the lower of the high and low frequency */ 

/* band magnitudes to a threshold level. If it is lower than the */ 

/* threshold the program updates the valid digit stack with a -1 */ 

/* and retrieves another frame. Otherwise, it continues and checks */ 

/* the twist. */ 

/* V 



rl=rl8+32 

a0= *rl - *rl + + 

rlO=_ttrlevl 

a2 = *rl 

al = *rl — - *rlO 

if (alt) goto _ttrlow 

r2=__twisthg 
_ttrhigh:a2 = *rl 

al = *rl++ - *rlO 

r2=_twistlw 

al = -*rl + a2 * *r2 

nop 

if (alt) goto notone 

nop 

if (alt) goto notone 

nop 

goto frqofset 

nop 
_ttrlow:al = -*rl + a2 * *r2 

if (alt) goto notone 

2*nop 

if (alt) goto notone 

nop 



/* loc . of max row & col mag's */ 
/* high mag. - low mag. */ 



/* lower value - 8e7 */ 

/* lower value - 8e7 */ 

/* low mag. * 2.5 - high mag. */ 

/* threshold check */ 

/* reverse twist test */ 

/* continue */ 



/* high mag. * 9.2 - low mag. 
/* threshold check */ 

/* normal twist test */ 



V 



/* */ 

/* Do largest frequencies stand out? */ 

/* This section compares the largest magnitude (detected "tone") in */ 

/* each frequency band with the other magnitudes in that band. If the */ 

/* "tone" is not significantly larger then the other magnitudes the */ 

/* program returns to begin and samples another frame. */ 

/* */ 



frqofset : r4=0 

rlO=_vldtonl 

rl=rl8+32 

r2=rl8 
freqoff 2:r3=2 
freqoffl:a2 = *r2 

al = *rl - *r2++ 

nop 

aO = *rl - a2 * *rlO 

nop 

if (aeq) goto freqoff 

nop 

if (age) goto freqoff 

nop 

goto notone 

nop 
freqoff: if (r3 — >= 0) goto freqoff 1 

nop 

rlO=_vldtonh 

if (r4 — >= 0) goto freqoff2 

rl=rl+4 



/* loc. of largest mag. in each band */ 
/* loc. of mag~2 of each tone */ 



/* check if mag's are the same */ 

/* highest mag. - other mags. * coeff */ 

/* if equal mag's get next mag */ 

/* if highest mag < other mag * coeff */ 
/* tone is not valid */ 

/* Therefore, update _ttdigit with -1 */ 
/* and get another frame of samples. */ 
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/* */ 

/* This section compares the total inband energy to the largest mag. */ 

/* in the low band, the largest mag. in the high band, and the sum of */ 

/* the two magnitudes. */ 

/* */ 



/* a2=largest mag in low band */ 

/* al=largest mag in high band */ 

/* a3=low band mag-energy*lowband coef */ 

/* a3=highband mag-energy*highband coef*/ 

/* a3=sum of mags-energy*sum coef */ 



energ:rl=rl8+32 
a2=*rl++ 
al=*rl 
r3=_nrgcoef 

a3 = a2 - *rl2 * *r3++ 
a2=al+a2 

a3 = al - *rl2 * *r3++ 
a3 = a2 - *rl2++ * *r3 
if (alt) goto notone 
nop 

if (alt) goto notone 
nop 

if (alt) goto notone 
nop 

/* */ 

/* Is the harmonic ratio > threshold? */ 

/* This computes the magnitudes squared of the "tones" second */ 

/* harmonics and compares it with the "tone". If the magnitudes are */ 

* within a certain ratio, the program returns to begin and samples */ 



/* another frame. */ 

/* */ 



r3=rl8+32 
rl = 
r5 = 

rlO=_iirincr 
r7=_hrmcoef 
_2ndharm:r6=*rl0++ 
r4=rl8 
r8=r6*2 
r8=r8*2 
r2=r8+32 
r2=r2+r5 
r8=r2*2 
r2=r2+_iirco 
r4=r4+r8 
a2= *r4 * *r4++ 
a l=* r 4 — * * r 4 
aO = al + *r4 * *r4 
aO = aO - a2 * *r2 
2*nop 

a3 = *r3 + + - aO * *r7++ 
3*nop 

if (ale) goto notone 
r5 = 16 

if (rl — >= 0) goto _2ndharm 
r9= iirincr 



/* loc . of largest mag. in each band */ 



/* compute the location of the 2nd harm */ 
/* coefficients and state variables */ 



/* loc. of the 2nd harmonic coeff's */ 

/* Sn * Sn-1 */ 
/* Sn-1~2 */ 
/* Sn~2 + Sn-1~2 */ 
/* |Y|~2 */ 

/* 1st harm - 2nd harm * 2nd harm coeff */ 



/* used in decode digit section */ 
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/* */ 

/* Decode digits */ 

/* pdr's lower byte contains the valid decoded digit in its lower 4 */ 

/* bits and the receiver number in its upper 4 bits. */ 

/* */ 

rl=*r9++ /* location of _iirincr */ 

r7=rl6*2 

r2=rl*2 

r2=r2*2 

r3=*r9 

r2=r2*2 

r3=r3*2 

r2=r2+r3 

r2=r2+_romtbl /* loc . of decoded digit */ 

rl0=*r2 /* present digit D(k) */ 

r4=rl9+2 

r5=*r4 /* 2nd to last digit D(k-2) */ 

r7=r7*2 

r5=r5-rl0 

if (eq) goto tonel /* D(k)=D(k-2) ? */ 

r4=rl9 

r5=*r4 /* last digit decoded D(k-l) */ 

r7=r7*2 

r5=r5-rl0 

if (ne) goto tonel /* D(k)=D(k-l) ? */ 

r8=r7*2 

r8=r8~rl0 

pdr=r8 

goto tonel 

nop 
notone:rlO=-l /* set digit to -1 and update stack */ 

tonel : rl=rl9 

r3=*rl 

*rl++ = rlO /* update D(k-l)=D(k) */ 

*rl = r3 /* update D( k-2 ) =D( k-1 ) */ 

rl9=rl9+4 

rl6=rl6+l 

if (rll — >= 0) goto timing /* process timing for next receiver */ 

rl8=rl8+128 

goto begin 

rl5=12 
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/* */ 

/* PROGRAM VARIABLES */ 

/* */ 

/* First order harmonic coefficients */ 

_iirco: float 1.70773781, 1.645281, 1.568687, 1.4782 
float 1.1641, .99637, .7986184, .5685327 

/* Second order harmonic coefficients */ 

float .916368, .70695, .46077885, .1851 

float -.644862, -1.0072464, -1.36221, -1.676771 

_iirzero: float 0.0 
_ttrlevl: float 4e7 
_twisthg: float 10.5 
_twistlw: float 7.4 
_vldtonl: float 5.0 
_vldtonh: float 5.0 
_nrgcoef: float 7.0, 5.6, 36.0 
_hrmcoef: float 1.6, 60.0 
_romtbl: int 1, 2, 3, 10 

int 4, 5, 6, 11 

int 7, 8, 9, 12 

int 14, 0, 15, 13 

/* */ 

/* */ 

/* STATE VARIABLES */ 

/* */ 

inbuf: 15*float 0.0 
ilimit: float 0.0 
dummy: 100*float 0. 
bufio: 16*float 0.0 
_iirincr : 2*int 
energy: 16*float 0.0 
_ttdigit:32*int 
.rsect ".hi_rai" 
_iirstvr:512*float 0.0 /* state variables */ 

/* */ 
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